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1 Introduction 



The physics of high-density QCD has become an increasingly active subject 
of research, both from experimental and theoretical points of view. In par- 
ticular, the collider facilities such as the BNL Relativistic Heavy Ion Collider 
(RHIC), and CERN Large Hadron Collider (LHC) will be able to probe new 
regimes of dense quark matter at very small Bjorken x or /and at large A, 
with rather different dynamical properties. Basically, new phenomena asso- 
ciated with an ultradense environment that may be created in the central 
collision region of these reactions are expected [|TJ. 

From the analysis of nucleus-nucleus collisions for RHIC energies and be- 
yond, we have that the perturbative QCD processes should determine the 
initial conditions, with most of the entropy and transverse energy presum- 
ably produced already during very early times (within the first 2 fm after 
the nuclear contact) by frequent, mostly inelastic, semihard gluonic collisions 
involving typical momentum transfers of only a few GeV 0. In particular, 
at early times, r 1/pr < 1/po ~ 0.1 fm for p « 2 GeV, semihard produc- 
tion of minijets will set the stage for further evolution of the system. The 
calculation of this process is based on the jet cross section for p T > p , with 
the parton densities evaluated at scale pt, with x values at central rapidities 
as low as x ~ O(10~ 4 ) in Pb + Pb collisions at LHC (5.5 TeV/nucleon). At 
the lower RHIC energies, x ~ O(l0~ 2 ) at central rapidities, and at higher 
rapidities the x values probed can be even smaller. Thus the small x behavior 
of the parton densities strongly influences the initial conditions of the minijet 
system. 

While the deep inelastic scattering data from HERA continues to refine 
the parton densities at small x, uncertainties in the distributions still exist, 
mainly associated to the high parton density effects present in this kinemat- 
ical region. Such effects would be present in nucleus-nucleus (AA) collisions 
at collider energies, modifying the perturbative QCD predictions for global 
observables, such as particle multiplicities and transverse energy production, 
as well as minijet production, heavy quarks, their bound states, and dilepton 
production |J. Consequently, the high density effects are one of the major 
theoretical issues in modeling the QCD processes in nuclear collisions. 

In this paper we analyze the high density effects in the behavior of the 
nuclear gluon distribution xGa and nuclear structure function F£ at small x 
and a large perturbative scale Q 2 . Our study is motivated by the perspective 
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that in a near future an experimental investigation of these effects at small 
x and Q 2 > 1 GeV 2 using eA scattering could occur in RHIC, as well as at 
DESY Hadron Electron Ring Accelerator (HERA). Furthermore, our goal is 
to improve the description proposed by Eskola, Kolhinen and Salgado (EKS) 
taking into account the high parton densities effects in the parton evolution 
at small x. 

In recent years several experiments have been dedicated to high precision 
measurements of deep inelastic lepton scattering (DIS) off nuclei. Experi- 
ments at CERN and Fermilab focus especially on the region of small values 
of the Bjorken variable x = Q 2 /2Mi > , where Q 2 = —q 2 is the squared four- 
momentum transfer, v the energy transfer and M the nucleon mass. The 
data [|], |5|, taken over a wide kinematic range, have shown that the pro- 
ton and neutron structure functions are modified by a nuclear environment. 
The modifications depends on the parton momentum fraction: for momen- 
tum fractions x < 0.1 and 0.3 < x < 0.7, a depletion is observed in the 
nuclear structure functions. The low x (shadowing region) and the larger x 
(EMC region) are bridged by an enhancement known as antishadowing for 
0.1 < x < 0.3. We refer to the entire phenomena as the nuclear shadowing 
effect. 

The theoretical understanding of F£ in the full kinematic region has pro- 
gressed in recent years, with several models which describe the experimental 
data with quite success ||. Here we will restrict ourselves to the descriptions 
which use the DGLAP evolution equations |7j] to describe the behavior of 
the nuclear parton distributions. Recently, Eskola, Kolhinen and Salgado 
[§], following Ref. f9j, have shown that the experimental results (§] present- 
ing nuclear shadowing effects can be described using the DGLAP evolution 
equations with adjusted initial parton distributions. The basic idea of this 
framework is the same as in the global analyzes of parton distributions in the 
free proton: they determine the nuclear parton densities at a wide range of 
x and Q 2 > Qq = 2.25 GeV 2 through their perturbative DGLAP evolution 
by using the available experimental data from I A DIS and Drell-Yan (DY) 
measurements in pA collisions as constraint. EKS have expressed the results 
in terms of the nuclear ratios Rf(x, Q 2 ) for each parton flavor / in a nucleus 
with A nucleons (A > 2), at 10~ 6 < x < 1 and 2.25 GeV 2 < Q 2 < WGeV 2 . 
The results of EKS seems to show that, in the kinematic region of the present 
data, the high density dynamical effects are small enough to be described by 
the DGLAP evolution equation with a suitable set of nonperturbative initial 
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conditions. The main shortcoming of use the EKS parameterizations in the 
calculations of eA, pA or AA processes is associated to the small x predic- 
tions, where the solution of the DGLAP equations reduces to the well known 
limit of the double logarithm approximation (DLA). This limit is character- 
ized by a strong growth of the gluon distribution, with a similar behavior 
for the F 2 structure function, which implies a high parton density in this 
kinematic region. However, when the density of gluons and quarks becomes 
very high the physical processes of interaction and recombination of partons, 
not considered in the DGLAP equations, become important in the parton 
cascade and these effects should be expressed in a new evolution equation. 
Therefore, the EKS description should be improved to include the high par- 
ton density effects when smaller values of x are considered. This is the case, 
for instance, in the calculation of the minijet cross section at LHC. 

At this moment, there are many approaches in the literature that propose 
distinct evolution equations for the description of the gluon distribution in 



high density limit M, 11[ |12], [13| . In general these evolution equations 



resum powers of the function k(x, Q 2 ) = ^\qi. A X9 ^§ ~ ; which represents the 
probability of gluon-gluon interaction inside the parton cascade. Moreover, 
these equations match (a) the DLA limit of the DGLAP evolution equation 
in the limit of low parton densities {k — > 0); (b) the GLR equation and the 
Glauber-Mueller formula as first terms of the high density effects. The main 
differences between these approaches are present in the limit of very large 
densities, where all powers of k should be resumed. Although the complete 
demonstration of the equivalence between these formulations in the region 
of large k is still an open question, some steps in this direction were given 
recently |T3|, |T5J. Here we will consider the Glauber- Mueller approach for 
the high density effects, which is a common limit of the current high density 
approaches in the kinematic region which we are interested. Thus we intend 
to obtain no model dependent predictions. 

The outline of this paper is the following. In next section we present 
a brief review of the Glauber-Mueller approach for the nuclear structure 
function (For details see In Sec. [|we analyze the EKS parameterization 
and present a procedure to improve this paramerization for the small x region, 
where the high density effects cannot be disregarded. Moreover, we present 
our results for the ratios xGa/xGn and F 2 A / ' and verify that at small x 
they are strongly modified in comparison with the EKS predictions. Finally, 
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in Sec. ^ we discuss the implications of our results for the heavy ion collisions 
in RHIC and LHC and present our conclusions. 



2 The High Density Effects in DIS 

The deep inelastic scattering off a nucleus is usually interpreted in a frame 
where the nucleus is going very fast. In this case the nuclear shadowing 
is a result of an overlap in the longitudinal direction of the parton clouds 



originated from different bound nucleons |]I0 |. Thus low x partons from 
different nucleons overlap spatially creating much larger parton densities than 
in the free nucleon case. This leads to a large amplification of the nonlinear 
effects expected in QCD at small x. In the target rest frame, the electron- 
nucleus scattering can be visualized in terms of the propagation of a small 
qq pair in high density gluon fields through much larger distances than it 
is possible with free nucleons. In terms of Fock states we then view the eA 
scattering as follows |16|]: the electron emits a photon (|e >— > \e"f >) with 
E 1 = v and pf ~ Q 2 , after the photon splits into a qq (\ery >— > \eqq >) and 
typically travels a distance l c ~ l/m^x, referred as the 'coherence length', 
before interacting in the nucleus. For small x (large s, where yfs is the 
center-of-mass energy), the photon converts to a quark pair at a large distance 
before it interacts to the target. Consequently, the space-time picture of the 
DIS in the target rest frame can be viewed as the decay of the virtual photon 
at high energy (small x) into a quark-antiquark pair, which subsequently 
interacts with the target. In the small x region, where x <C the QQ 

pair crosses the target with fixed transverse distance r t between the quarks. 



Following Gribov Jl6||, we may write a double dispersion relation for the 



forward elastic amplitude A, related to the total cross section by the 
optical theorem [ImA = sa(s, Q 2 )}, and obtain for fixed s 

r dM 2 dM' 2 1 
a(s ' g2) =?/ W + Q 2 M' 2 + Q 2 p ^ M2 ' M ' 2) S Im A ^ M2 ' M ' 2) ' ( 

where M and M' are the invariant masses of the incoming and outgoing 
qq pair. If we assume that forward qq + A scattering does not change the 
momentum of the quarks then A q g + A is proportional to 5(M 2 — M' 2 ), and 
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flip becomes 

r dM 2 

Q 2 ) = E J (M2 + g2)2 P(^ M 2 )a q - q+A (s, M 2 ) , (2) 

where the spectral function p(s, M 2 ) is the density of gg states, which may be 
expressed in terms of the 7* — > gg matrix element ||17|| . Using that M 2 = (fc 2 + 



m 2 )/[z(l — z)], where fcf and z are the transverse and longitudinal momentum 
components of the quark with mass m q , we can express the integral over the 
mass M of the gg in terms of a two-dimensional integral over z and k t . Instead 
of kt is useful to work with the transverse coordinate r t (impact parameter 
representation), which is the variable Fourier conjugate to k t , resulting [JIB 



F 2 A (x,Q 2 ) = -^a(s,Q 2 ) 

47TQ; em 

Q 2 f J fd 2 r t 



dz [ —\^(z,r t )\ 2 a qq+A (z,r t ) , (3) 

J TV 



where 



n)\ 2 = ^ E e){[z 2 + (1 - ,) 2 ]e 2 Kl (er t ) 2 + m 2 K (er t ) 2 } . (4) 

The photon wave function ^f(z, r t ) is simply the Fourier transform of the ma- 
trix element for the transition 7* — > gg. Moreover, a em is the electromagnetic 
coupling constant, e 2 = z(l — z)Q 2 + m 2 , mi is the quark mass, n/ is the 
number of active flavors, e 2 is the square of the parton charge (in units of e), 
Kq^ are the modified Bessel functions and z is the fraction of the photon's 
light-cone momentum carried by one of the quarks of the pair. In the leading 
log(l/x) approximation we can neglect the change of z during the interaction 
and describe the cross section a q ^(z, 4:/r 2 ) as a function of the variable x. To 
estimate the high density effects we consider the Glauber multiple scattering 
theory |19| |, which was derived in QCD f20|| . In this framework the nuclear 



collision is analyzed as a succession of independent collisions of the probe 
with individual nucleons within the nucleus, which implies that 

F 2 A (x,Q 2 ) = (dz I ^\M,( z ,r t )\ 2 [ fh 2 [l- e ~^ +N ^ s ^} , (5) 



47ra em J J it J it 
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where bt is the impact parameter, S(bt) is the profile function and a qq+ is 
the dipole cross section off the nucleons inside the nucleus, which is propor- 
tional to the pair separation squared rf and the nucleon gluon distribution 
xg{x,l/rj). The expression fl5|) represents the Glauber-Mueller formula for 



the nuclear structure function (see |TT for details). The use of a Gaussian 



parameterization for the nucleon profile function S{pt) = ^-e r a , where 
Ra is the mean nuclear radius, simplifies the calculations. We obtain that 
the F2 structure function can be written as 



F 2 A (x, Q 2 ) = ^E A jf° + InMx, rl)) + El ( Kq (x, rf))} , (6) 

^ u,d,s Cp ^ t 

where n q = (2a s A/3R 2 ) 7rrf xGn(x, ^), C is the Euler constant, E\ is the 
exponential integral function and A the number of nucleons in a nucleus. 
This equation allows to estimate the high density corrections to the structure 
function in the DLA limit. Expanding the equation (§) for small K q , the first 
term (Born term) will correspond to the usual DGLAP equation in the small 
x region. 

The Glauber - Mueller formula has been used in a comprehensive phe- 
nomenological analyzes of the behavior of distinct observables in ep and eA 
processes. The results of these studies agree with the current ep HERA data 



^] and allow us to make some predictions which will be investigated 
in the near future in eA colliders 0, p6|l . The main shortcoming in the 
studies of eA processes is that the large x effects (the antishadowing and 
the EMC effect) were disregarded, which implies predictions only for the 
asymptotic behavior (large s/small x) of the observables. In next section we 
propose a procedure to obtain more realistic predictions in the full kinematic 
region. 



3 The improved EKS parton distributions 

As discussed in Sec. [1] the EKS parameterization, although describes the 
current experimental fixed target data quite well, is not a good approxima- 
tion for the small x region in the perturbative regime, where the DGLAP 
equations predicts a strong growth of the parton distributions (the DLA 
limit). In this limit the nuclear gluon distribution is given by xG a oc 



6 



exp y/ln(l/x) ln(Q 2 /(5o)> which is almost independent of the initial nonper- 
turbative input, i.e. of the adjusted parameters obtained in the EKS param- 
eterization. To improve the EKS description we propose the substitution of 
DLA limit of the DGLAP evolution, present in the EKS framework, by the 
Glauber-Mueller evolution. This procedure implies to calculate the nuclear 
structure function using the following expression 

F A (x,Q 2 ) = F A (x,Q 2 )[EKS] ~ (l/A)F A (x,Q 2 )[DLA] 

+ (l/A)F A (x,Q 2 )[eq. (6)], (7) 

with F 2 A (x, Q 2 )[EKS] = Rp 2 x F 2 N (x, Q 2 ), where Rp 2 is obtained in terms of 
a combination of nuclear parton ratios R A (x,Q 2 ) [[J. The nucleon structure 
function is given by F 2 N (x, Q 2 ) = J2u,d,s e q [xq(x, Q 2 ) +xq(x, Q 2 )} +F 2 (x, Q 2 ), 
where the charm component of the nucleon structure function is calculated 



considering the charm production via boson-gluon fusion mechanism |21| and 
the nucleon parton distributions are given by the GRV parameterization |H] . 
In this work we assume m c = 1.5 GeV. Moreover, F A (x, Q 2 )[DLA] is cal- 
culated using the parton distribution from DLA limit of DGLAP evolution. 
For practical purpose, this term is given by the Born term of Eq. (^|). 

The above procedure implies the inclusion of: (a) the full DGLAP evo- 
lution equation in all kinematic region; (b) the nonperturbative nuclear cor- 
rections in the nuclear parton distributions, that describes the experimental 
fixed target data, and (c) the high density effects present in the parton evo- 
lution at small x in the perturbative regime. Therefore, with the Eq. ([?]) we 
are able to describe the parton evolution in all kinematic region. 

A similar description can be used to estimate the nuclear gluon distribu- 
tion xG A . In this case we take 

xG A (x,Q 2 ) = xG A (x,Q 2 )[EKS] - (l/A)xG A (x,Q 2 )[DLA} + 

+ (l/A)xG A (x,Q 2 )[GM] (8) 

where xG A (x, Q 2 ) [EKS] = Rq x xG n (x,Q 2 ) is the EKS prediction, 
xG A (x, Q 2 ) [DLA] is the DGLAP (DLA) prediction for the nuclear gluon dis- 
tribution and GM represents the Glauber-Mueller nuclear gluon distribution 
given by [T| 

xG A (x,Q 2 )[GM] =ff^ ^{C + lniKGix'y^ + E^KGix'y,))} 

it z Jx x' ii Trry 
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where kq = (9/4)k q . The main difference between the high density effects 
in the quark and gluon densities stems from the much larger cross section 
a N = (9/4) cjy, i. e. = (9/4)/t g , which in turn leads to a much larger 
gluon shadowing. As before, the DLA limit is obtained from the Born term 
of the Glauber-Mueller evolution [Eq. (J|)]. 

Using the above procedure we can calculate the ratios R F2 = F 2 A / F 2 D and 
Rg = xGa/xGn, where F 2 D is the deuterium structure function. In order to 
include only perturbative contributions to the estimative of the high density 
effects, we take the initial scale at Q$ = 1 GeV 2 in the Glauber- Mueller 
expressions [Eqs. (^) and (§)] and calculate the ratios for A = 208. 

In Fig. |l] we present our predictions for the behavior of the ratio R F2 = 
F2/F2 (denoted as EKS MOD) for two values of virtualities (Q 2 = 2.25 and 
15 GeV 2 ). For comparison the EKS predictions are also shown. We verify 
that while for the region of large values of x (> 10~ 2 ) our predictions are 
almost identical to the EKS results, at small x the differences are large. A 
comment related to experimental data is in order. As x ~ Q 2 /s, where s is 
the squared CM energy, the data in the region of small x values are for small 
values of Q 2 (< 1 GeV 2 ), where the use of the perturbative QCD cannot be 
justified and the shadowing corrections are dominated by soft contributions. 
Therefore, our perturbative predictions cannot be compared with the current 
experimental fixed target data. In the perturbative regime Q 2 > 1.0 GeV 2 , 
where the data are associated with x values greater than 10~ 2 , our results 
are almost identical to the EKS predictions. Thus, the fixed target data does 
not allow us to discriminate between the predictions. 

The results shown in Fig. [I] predicts a sizeable modification in the quark 
distribution at low x (x < 10~ 3 ). Thus, we expect that the I A cross section 
will present a strong reduction when compared to the Ip in this kinematic re- 
gion. Also the Drell-Yan production in heavy ion collisions should be strongly 
modified when compared to the pp one at LHC energy (x ~ 10~ 4 ). At RHIC 
kinematic region, the high density effects does not seem to be strong enough 
to modify the Drell-Yan cross section in respect to the EKS prediction. Nev- 
ertheless, this subject deserves a further more detailed study. 

The high density effects are important already for the initial scale of the 
EKS parameterizations (Q 2 = 2.25 GeV 2 ) and increases with the virtuality. 
The remarkable feature in the improved parton distributions is the nonsat- 
uration of the ratio at small values x. In the EKS parameterization the 
saturation is included in the initial condition and this general behavior is not 
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modified by the evolution. This is a consequence of the DGLAP evolution 
equations, which reduces to the DLA limit at small x in the nuclear and nu- 
cleon case, keeping the ratio approximately constant. When the high density 
effects are considered the ratio is not constant since these effects increases 
at small values of x and are larger in the nuclear case. Therefore, the non- 
saturation of the ratio in the perturbative regime is a signature of the high 
density effects in the nuclear processes. 

In Fig. H] we present our predictions for the behavior of the ratio Rg(x, Q 2 ) = 
xGa/xGn (denoted as EKS MOD) for two values of virtualities (Q 2 = 2.25 
and 15GeV 2 ). For comparison the EKS predictions are also shown. We 
verify that while for the region of large values of x (> 10~ 2 ) our predictions 
are almost identical to the EKS results, at small x the differences are large. 
We can see from the both figures that the EKS assumption that Rq = R F2 
at small values of x is strongly modified by the high density effects and its 
Q 2 evolution, since Rq -C R F2 for Q 2 = 15GeV 2 . Therefore, we predict a 
large modification in the quarkonium production at LHC, as well as in its 
bound states (e.g. J/^, T, ...). Furthermore, as the high density effects 
are important already to Q 2 ~ 2.25 GeV 2 , also the minijet production will be 
modified at LHC. At RHIC kinematical region the high density effects does 
not significantly modify the gluon distribution and we expect a similar pre- 
diction of the quarkonium production from both models. Anyway, the high 
density effects should be considered in detail before to use the quarkonium 
production, as well as its bound states, as probes of the deconfined state of 
matter. Also this subject deserves a further detailed study. 

It is important to note that the EKS description use as a constraint the 
momentum sum rule. We verify that the improved distributions violates this 
sum rule at must in 5 % for the gluon distribution, which is small when com- 
pared with the experimental uncertainty on this distribution in the antishad- 
owing region. Moreover, the violation is small since the main contribution 
to the sum rule comes from the large x region, where the high density effects 
are negligible. 

4 Conclusions 

In this paper we have proposed an improvement of the EKS nuclear parton 
distributions in the small x region by the inclusion of the perturbative high 
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density effects. We predict the nonsaturation of the ratios R Fi and Ro when 
these effects are considered and that R F2 <C Rg in the small x limit. Such 
results could be tested in the future eA colliders. Furthermore, our results 
demonstrate that the high density effects are very important mainly at the 
LHC kinematic region, where strong modifications in the Drell- Yan and 
quarkonium production are expected. As the small x behavior of the parton 
densities strongly determines the initial conditions of the minijet system in 
nucleus-nucleus collisions, our results show that the high density effects can- 
not be disregarded in the calculations of the observables and signatures of a 
Quark-Gluon Plasma. 
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Figure Captions 



Fig. [l]: Comparison between the predictions of the original EKS and EKS 
modified (denoted EKS MOD. in plot) for the ratio R F2 (x, Q 2 ) = Ff/Fj> as 
function of the variable x at different values of Q 2 . See text. 

Fig. § Comparison between the predictions of the original EKS and EKS 
modified (denoted EKS MOD. in plot) for the ratio Rg(x,Q 2 ) = xGa/xGn 
as function of the variable x at different values of Q 2 . See text. 
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Figure 2: 



15 



